Draft February 2, 2008; for submission to ApJL 



The Spitzer Gould Belt Survey of Large Nearby Interstellar 
Clouds: Discovery of a Dense Embedded Cluster in the 
, Serpens- Aquila Rift 

O 

^ : R. A. Gutermuthi, T. L. Bourke\ L. E. Allen\ P. C. Myers\ S. T. Megeath^ B. C. 

Q ■ Matthews^, J. K. J0rgensen^, J. Di Francesco^, D. Ward-Thompson^, T. L. Huard^, T. Y. 

^ '. Brooke^ M. M. Dunham^ L. A. Cieza^ P. M. Harvey^ N. L. Chapman* 



ABSTRACT 

Q^- We report the discovery of a nearby, embedded cluster of young stellar objects, 

Q . associated filamentary infrared dark cloud, and 4.5 /im shock emission knots from 

^ ■ outflows detected in Spitzer/IRAC mid-infrared imaging of the Serpens-Aquila 

ci . Rift obtained as part of the Spitzer Gould Belt Legacy Survey. We also present 

radial velocity measurements of the region from molecular line observations ob- 
^ . tained with the Submillimeter Array (SMA) that suggest the cluster is co-moving 

2 . with the Serpens Main embedded cluster 3° to the north. We therefore assign 

, it the same distance, 260 pc. The core of the new cluster, which we call Ser- 

pens South, is composed of an unusually large fraction of protostars (77%) at 



^ ' high mean surface density (>430 pc~^) and short median nearest neighbor Spae- 
th ! ing (3700 AU). We perform basic cluster structure characterization using nearest 
^ ' neighbor surface density mapping of the YSOs and compare our findings to other 
- , , , known clusters with equivalent analyses available in the literature. 
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Subject headings: ISM: individual (IRAS 18275-0203, IRAS 18274-0205) - stars: 
formation - stars: low-mass 



Introduction 



The aim of the Spitzer Gould Belt (GB; Allen et al. in prep.) and c2d (lEvans et al. 



20031 ) Legacy Surveys is to build a complete mid-infrared record of star formation in all 
of the large molecular clouds within 500 pc, largely dominated by the Gould Belt itself. 
This is one of several surveys of the Gould Belt bei ng performed at many wavelengths from 
the near-IR (Two- Micron All Sky Survey; 2MASS ISkrutskie et al.ll2006l ) through the far- 
infrared {Herschel; lAndre fc Saracend l2005r) to the submillimeter (Submillimeter Common 
User Bolometer Array 2; SCUBA-2; IWard- Thompson et al.l 120071 ). In combination, the re- 
sulting database offers a nearly complete and unbiased view of all nearby star formation in 
large molecular clouds. This penetrating examination of regions that have often received 
little attention in the literature is likely to reveal new regions of star formation. 

The Spitzer Infrared Array Camera (IRAC) imaging presented here reveals a dense 
cluster of sources with excess infrared emission and associat ed 4.5 //m-bright structured 
nebulosity (typically shocked H2 emission from outflows, cf. ISmith et al.l l2006l ) within a 
filamentary infrared dark cloud seen in absorption against bright diffuse 5.8 a nd 8.0 fim 
emission. This emission is most lik ely PAH-feature em ission ( iMorris et al.ll2004l ) excited by 
the young, high mass stars in W40 (ISmith et al.lll985[ ). 20' to the east. We found structured 
Ks nebulosity at the site in 2MASS Atlas images, but few sources. The region is flanked by 
two Infrared Astronomy Satellite (IRAS) sources (IRAS 18275-0203 and IRAS 18274-0205), 
and Midcourse Space Experiment (MSX) coverage of this region reveals only a few additional 
point sources. No deep X-ray images of this region are available. 



2. Observations 



We obtained 3.6, 4. 5, 5.8, and 8.0 /x m imaging of 5.1 square degrees of th e Serpens- 
Aquila Rift with IRAC (IFazio et al.ll2004l ) onboard Spitzer ( iWerner et al.l 120041 ) as part of 
the GB Legacy Survey (PID: 30574) on 27 October 2006. In this letter, we focus on the 
14' X 10' field of view centered on Serpens South (R.A., Dec. (2000) = 18:30:03, -02:01:58.2). 
High Dynamic Range mode was used, resulting in 0.4 and 10.4 second integrations for each of 
four dithered images at each mosaic position. Mosaics were constructed at 1'.'2 per pixel with 
Basic Calibrated Data (BCD) products from the S]9zteer Science Center's data pipeline vl4.4. 
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All artifact treatment, mosaic construction, point source detection and aperture photometry, 
and photometric classification of sources pr esented here has been p erformed using custom 
Interactive Data Language (IDL) routines (iGutermuth et al.l 120071 ). We adopted aperture 
radii and inner and outer background annulus ra dii of 2'! 4, 2M , and 7^/2, respectively. Point 
source flux calibration is discussed in detail in iReach et al.l (120051 ). Field averaged 90% 
differential completeness magnitude limits are 14.5 , 14.2, 13.2, and 12.0 f or the 14' x 10' 
submosaics at 3.6, 4.5, 5.8, and 8.0 fim, respectively (IGutermuth et al.ll2005r). The se nsitivity 
of these data are equivalent to the IRAC data presented by IWinston et al.l (120071 ) for the 
Serpens Main cluster, which will be discussed in Section |H IRAC photometry catalogs are 
bandmerged together with a radial matching tolerance of 1", and then merged with the 
2MASS Point Source Catalog (PSC) at a tolerance of 1". 



2.1. Submillimeter Array 

Observations with the SM7^0 (jHo et al. 20041 ) were performed on 30 January 2007. The 
array was tuned to cover the frequency range 219.4-221.4 GHz (LSB) and 229.4-231.4 GHz 
(USB), which includes the lines of ^^CO, ^^CO, and C^^O 2-1, with a resolution of ~0.5 
km s~^ /channel. Three positions were observed for 30 minutes each; IRAS 18275-0203 
(R.A.,Dec. (2000) = 18:30:05.8, -02:01:45) and a nearby embedded source (R.A.,Dec. (2000) 
= 18:30:01.3, -02:01:48, labelled "Embedded" in Figure [ID, both in Serp ens South, and a 



positi on near to the CO 1-0 peak of W40, as determined from the maps of iGrabelsky et al. 



(119871 ). at R.A.,Dec. (2000) = 18:31:20.7, -02:01:07. Ganymede was observed for passband 
calibration. No flux or complex gain calibration was performed. More detailed SMA results 
will be discussed in a future paper. 



3. Results 



3.1. Distance 



Serpens South appears to be associated with a filamentary cloud seen in absorption 
against PAH emission from W40 (Fig . [J). Thus the distance to W40 (550 ± 150 pc 



Radhakrishnan et al. 1972; Smith et al. 



19851 ) is an upper limit to the distance to Serpens 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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South. The Serpens Main cluster has local standard of rest (LSR) velocit ies in the range 
6-11 kms~^, when observed in the isotopologues of CO (IWhite et al.lll995l ). with Gaussian 
FWHM widths of ~ 2 kms~^. Thus, its molecular line emission can be separated from that 
associated with W40, which occu rs in the LSR range of 2-6 kms~^ (IZeilik fc Ladal Il978 



Vallee et al.lll992l : IZhu et al.ll2006l ). and is spatially compact. The SMA data clearly show 



that the velocities of IRAS 18275-0203 and the nearby embedded source are the same as that 
of the Serpens Main cluster, and no t of the W40 cluster. We therefore assign the distance 
of 260 ± 37 pc (IStraizys et al.lll996l ) to S erpens South. This i s consistent with more recent 
results for the entire Serpens- Aquila Rift (IStraizys et al.ll2003l ). 



3.2. Identification, Membership, and Spatial Distribution of YSOs 



We identify likely YSOs from their excess infrared emission over that expect ed from a 

typic al photosphere and high infrared flux compared to extragalactic contaminants (IGutermuth et al 
20071 ). These candidates are classified as Class I (protostar) or Class II (pre-main sequence 
star with disk) YSO^by the their value of ajRAc^ the linear least squares f it to their spectra l 



energy distribution (log \F\ vs. log A) through the four IRAC bandpasses (ILada et al 



2006|). 



In this way, we ha ve identified 54 Class I and flat spectrum "protostars" {aiRAc > —0.3 
Greene et al1ll994h and 37 Class II (-1.6 < amAc < -0.3) YSOs, a total of 91 (59% pro- 
tostars), within the 14' x 10' (1.1 pc x 0.8 pc) field of view presented in Fig. [H Half of the 
sources in this field of view are located in a region that stands out as both relatively high 
in surface density and dominated by sources that are protostellar. Of the YSOs identified, 
37 Class I and 11 Class II sources are located within a 2'.5 (0.2 pc) radius circle centered on 
R.A., Dec. (2000) = 18:30:03, -02:01:58.2, making up the core of the Serpens South Cluster 
(Fig. d]). In this core, 77% of the YSOs are protostellar, with a mean density of 430 pc^^. 
This dense grouping appears elongated in an alignment similar to the dust filament seen in 
absorption (150° east of north). The asymmetric distribution suggests that the mean surface 
density of 430 pc~^ in the core is likely a significant underestimate. To demonstrate this, 
we note that the median projected distance between nearest neighbor YSOs here is 13'.'2, or 
3700 AU. Sources uniformly spaced at a density of 430 pc~^ are 21'.'4 apart, or 5600 AU. 



^ Class III (diskless) YSOs cannot be distinguished fr om field stars in this w ay (jMegeath et al. 
can be efi^ectively identified with X-ray observations (cf. iFeigelson et al.ll2007l ). 



200J), but 
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3.3. YSO Surface Density Mapping 



To characterize the two dimensional structure of the Serpens South cluster, we have 
constructed a nearest neighbor surface density map of all t he ^'pifeer-identified YSO s (Fig. [2]). 
The method used to generate these maps is documented in Gutermuth et al. ( 2005l N. though 
we have used n = 11 nearest neighbors her e to sample the surface densities at higher fidelity 
(33% uncertainty; ICasertano fc Hutlll985l ). The elevated density is evident in the center of 
the map, with a peak density of 1600 pc^^. The apparent boundary contour for the dense 
cluster core is 590 pc~^, a more reasonable mean surface density than was measured over the 
poorly matched circular area used in Section \^72\ above. Lower density star formation along 
the dust filaments is also apparent at densities from 50 to 120 pc~^. The overall structure 
is elongated, refiecting the same orientation and extent as the filamentary cloud structures 
seen in absorption (Fig. [1]). 



4. Putting Serpens South in Context 



The Serpens South core's high mean surface density (>430 pc~^), exceptionally large 
fraction of Class I sources (77%), and relatively large membership (48, with 43 more along 
the 0.5 pc cloud filaments) suggest both a very recent onset o f star formation (within the 
typical lifetime of the protostellar phase, 2 ± 1 x 10^ yr, cf. iKenyon &: HartmannI Il995l ) 
and high star formation rate (~90 Mq Myr~^ in the core, assuming the above protostellar 
phase lifetime and 0.5 M f^ per source) . Visually, the region bea rs some resemblance to 



the Serpens Main cluster (IHarvey et al.l l2006l : IWinston et al.l 120071 ) that we have shown lie 



at similar distances. In Fig. [3l we present 6.5' x 6.5' (0.5 pc x 0.5 pc) 8 fim images of 
the dense cores of both the Serpens South and Serpens M ain clusters, overplot ted with the 
positions of Spitzer identified Class I and Class II sources (IWinston et al.ll2007l ). Over these 
fields of view, the two regions have nearly identical numbers of YSOs with excess infrared 
emission and similar protostar fractions (55 at 71% and 54 at 57% in Serpens South and 
Serpens Main, respectively). However, Serpens South is more concentrated and confined 
to its filament axis than Serpens Main, as evidenced by median nearest neighbor distances 
between YSOs of 3700 AU and 4800 AU, respectively. These results suggest that Serpens 
South is at least as young and as prolific a star forming region as Serpens Main, though 
sampling statistics prevent us from making any claims beyond similarity. At larger scales. 



•^In summary, from each position in a imiform grid, we measure the distance j) that defines a circle 
that contains the nearest n sources. From this radial distance, a surface density is directly computed as 



= ^f^^ (ICasertano fc Hutlll985D . 
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both Serpens Main and Serpens South are active star forming sites embedded in larger, less 
dense, less extinguished distributions of YSOs. Given the high likelihood that both clusters 
are associated with the Serpens-Aquila Rift and their structural and evolutionary similarity 
over such a large projected distance from each other (197' or 15 pc), we speculate that the 
initial conditions of this region play a strong role in defining these qualities. 

Serpens South joins a growing number of young clusters that are protostar rich and 
asymmetrically structured. The more typical asymmetric clusters (e.g. NGC 1333, IRAS 20050+2720, 
GGD 12-15) have evolved to the point that star formation is active throughout most sites 
of dense gas in their natal clouds. This is supported by their large numbers of protostars 
(active current star formation), their larger numbers of stars with disks (star formation in 



( Gutermuth 


2005; 


Allen et al. 


2007) 



a particularly hi gh protostar fraction th ough (59%; NGC 1333 is 29% protostars, excluding 



transition disks; iGutermuth et al.l 120071 ) . and it is part of a dense, dusty cloud filament that 
is currently forming stars at a rather low surface density compared to the core. If the rest of 
the filament continues forming stars within the next few 10^ yr, we speculate that this region 
could evolve into an asymmetric cluster of similar number of sources and protostar fraction 
to the aforementioned asymmetric examples. Therefore, the Serpens South cluster could be 
described as a protocluster, a region with a sufficiently high star formation rate to produce 
a large number of stars before gas dispersal processes (outflows, radiation, etc.) from the 
growing cluster membership terminate the process. Given its close proximity, follow-up ob- 
servations (e.g. to characterize the gas distribution and kinematics) and detailed analysis 
should offer a new and unique view of the earliest stages of clustered star formation. 
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Right Mcension 

Fig. 1. — At left, a color-composite image of the IRAC mosaics of Serpens South; 3.6, 
4.5, and 8.0 fim images are mapped to blue, green, and red, respectively. At right, the 
spatial distribution of YSOs overlaid on the grayscale 8.0 /im image. Red circles are Class I 
protostars, and green diamonds are Class II stars with disks. The white circle marks the 
dense core of the cluster. 
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Right Aso^r^^SoHi 



Fig. 2. — A n = 11 nearest neighbor surface density map for the Serpens South YSOs. 
Contour levels mark stellar surface densities of 1 cr (33%) below successive contours (1300, 
890, & 590 pc~^ in white and 400, 260, 180, 78, & 52 pc~^ in black), and the linear inverse 
grayscale sets 1000 pc~^ as black and pc~^ as white. The field of view is the same as for 

Fig.m 
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Se^'pens Worn Core 




Fig. 3. — Spitzer/IRAC 8 yum im ages of the 0.5 pc x .5 pc core regions of the Serpens South 



(this work) and Serpens Main (jWinston et al 



20071 ) clusters. Overlaid are the positions 



of the Class I sources (circles) and the Class II sources (diamond s) for each region. We 
have included the fiat spectrum sources from IWinston et al.l (120071 ) as part of the Class I 
population for the purposes of this plot. 



